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Abstract: The aim of this paper is to study the existence and approximation of periodic solutions for non-linear 
systems of integral equations, by using the numerical- analytic method which were introduced by Samoilenko[ 
10, 11]. The study of such nonlinear integral equations is more general and leads us to improve and extend the 
results of Butris [2]. 
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I. Introduction 

Integral equation has been arisen in many mathematical and engineering field, so that solving this kind 

of problems are more efficient and useful in many research branches. Analytical solution of this kind of 

equation is not accessible in general form of equation and we can only get an exact solution only in special 

cases. But in industrial problems we have not spatial cases so that we try to solve this kind of equations 

numerically in general format. Many numerical schemes are employed to give an approximate solution with 

sufficient accuracy [3,4,6, ,12,13,14,15]. Many author create and develop numerical-analytic methods [1, 2,5, 

7,8,9] and schemes to investigate periodic solution of integral equations describing many applications in 

mathematical and engineering field. 

Consider the following system of non-linear integral equations which has the form: 
t 

x(t,x 0 ,y 0 ) = F 0 (t) + j /i(s,x(s,x 0 .yo).y(s.*o.yo). 
o 

s 

, j G 1 (s,T)g 1 (T l x(T l XQ,yQ) l y(T,XQ,yQy)dT, 

-co 

' C 9i(j,x(T,x 0 ,y 0 ),y(j,x 0 ,y 0 ))dT)ds •••( Ij ) 

t 

y(t,x Q ,y 0 ) = G 0 (t) + j f 2 (s,x(s,x 0 ,y 0 ),y(s,x 0 ,y 0 ), 
o 

s 

, j G 2 (s,T)g 2 (T,x(T,x 0 ,y 0 ),y(j,x 0 ,y 0 ))dT, 

— CO 

b(s) 

, j g 2 (T,x(T,x 0 ,y 0 ),y(T,x 0 ,y 0 ))dT)ds -(I 2 ) 

a(s) 

where x 6 D c R n ,D is closed and bounded domain subset of Euclidean space R n . 

Let the vector functions 
A (t, x, y, z, w) = (/ n (t, x, y, z, w), f 12 (t, x, y, z, w), ... , f ln (t, x, y, z, w)), 
f 2 (t, x, y, u, v) = (f 21 (t, x, y, u, v), f 22 [t, x, y, u,v),..., f 2n (t, x, y, u, v)), 
gi(t,x,y) = (g n (t,x,y),g 12 (t,x,y),...,g ln (t,x,y)), 
g 2 (t,x,y) = (g 21 (t,x,y),g 22 (t,x,y),...,g 2n (t,x,y)), 

F 0 (t) = (F 01 (t),F 02 (t) F 0n (t)), 

and 

G 0 (t) = (G 01 (t),G 02 (t),-,G 0n (t)) 

are defined and continuous on the domains: 

(t,x,y,z,w) ER 1 xD xD 1 xD 2 = (-oo, oo) x D x D l x D 2 x D 3 ) 
(t, x,y,u,v) E R 1 x D x D 1 x D 2 — (-oo, oo) x D X D 1 X D 2 X D 3 ) 
and periodic in t of period T, Also a(t) and b(t) are continuous and periodic in t of period T. 

Suppose that the functions f^t.x.y, z, w), f 2 (t,x,y,u, v), gxit.x.y) 
and g 2 (t, x, y) satisfies the following inequalities: 
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\\f 1 (t,x,y,z,w)\\ < M x , \\f 2 {t,x,y,u, v)\\ <M 2 , M 1 ,M 2 >0] 

Il5i(t.^y)ll <Af 1( (t.*.y) II <Wi , n 1 ,n 2 >o}"' { ■ ) 

Wfdt.xuyuzuwy) - f 1 (t,x 2 ,y 2 ,z 2 ,w 2 )\\ < K t \\x t -x 2 \\ + tf 2 llyi-y2ll+ 

+K 3 \\ Zl -z 2 \\ +K, || Wl -w 2 1| -(1.3) 
||/ 2 (t^y^Ui,!?!) -/ 2 (t,x 2 ,y 2 ,u 2 ,i7 2 )|| < Lilki -x 2 || + L 2 ||yi-y 2 ||+ 

+L 3 ll"i-"2ll +i 4 lbi-V2ll . -(1-4) 

llsi(t.*i.yi) -5i(t.^2.y 2 )ll < ffilki -* 2 II + ff 2 llyi -y 2 ll . -(1.5) 

\\9i{t,x Y , yx ) - g 2 {t,x 2 , y2 )\\ <H x \,x x -x 2 \\ +H 2 \\ yi -y 2 || , -(1.6) 

for all tSR 1 , x,x 1} x 2 £ D ,y , yi ,y 2 £ D 1 ,z,z 1 ,z 2 ,u,u 1 ,u 2 £D 2 , and w, w x w 2 , v, v x , v 2 ££> 3 , where M x = 
(M n ,M 12 M ln ), 

M 2 = (M 21 ,M 22 M 2n ), N, = (N n ,N 12 JV ln ) and 

N 2 = (JV 21 ,JV 22 , ...,N 2n ) are positive constant vectors and K t = (K Uj ), K 2 = (K 2ij ), K 3 = (K 3ij ), K 4 = (K 4ij ), 
L 1 = {L Uj \L 2 = (L 2ij ), 

L 3 = (L 3i . ), L 4 = (L 4y ), Ri = (%•). Ri = (.R 2iJ ), Hi = (%), and 
H 2 = (H 2i j ) are positive constant matrices. 

Also G 1 (t, s) and G 2 (t, s) are (nXn) continuous positive matrix and 
periodic in t, s of period Tin the domain R 1 X R 1 and satisfy the following conditions: 
IIGiCt,^)!! < y e -^iCt- s ) , y , A 1>0 -(1.7) 
and 

||G 2 (t,5)||<5e- 1 2( t - s ) , 5, A 2 >0 -(1.8) 
where -oo < 0 < s < t < T < oo , i,j = 1,2, — ,n, 
andh = max t£[0 T] \b(t)-a(t)\ , \\ . \\ = max t€[0J] \ . \ 



Now define a non-empty sets as follows: 



Df = D — —M 1 , 



D 



1/2 



-M ■ 
2 

T/Y 



(1.9) 



D 2fl =D 2 --(P)(R 1 M 1 +R 2 M 2 ) , 
T 

%! =D 3 --h(R 1 M 1 + R 2 M 2 ) , 
D 2f2 =D 2 - T - (-) Qi 1 M 1 + H 2 M 2 ) , 



D 3h =D 3 --h{H 1 M 1 +H 2 M 2 ) 

Forever, we suppose that the greatest eigen-value q max of the following 

matrix: 

Qo 

is less than unity, i.e. 



^ [L, + H 1 (L 3 |- + hi,)] T - [L 2 + tf 2 (|-L 3 + hi,)} 



\ 



<lmax (<2o) 

■Y 



<Pi + VV + 4 (^2 -<Ps) 



< 1 



(1.10) 



(1.11) 



where ^ = \ \K X + fl^tf, + hK,)] +\[L 2 + H 2 (^L 3 + ftL 4 )], 
cp 2 = (£[/f 2 + R 2 (^K 3 + WQ)])(£[£i + #1(^3 + ftL 4 )]) and 
<p 3 = ^[K t + ffiC^/C, + Wf 4 )])(j[L 2 + H 2 (^L 3 + «. 4 )]). 

By using lemma 3.1 [10], we can state and prove the following lemma: 
Lemma 1.1. Let /i (t, x, y, z, w) and f 2 (t, x, y, u, v) be a vector functions which are defined and continuous in 
the interval [0, T], then the inequality: 
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/||Li(t,x 0 ,yo)ll\ < (M 1 a(t)\ 

\\\L 2 (t,x 0 ,y 0 )\\)-{M 2 a(t)) "' U1ZJ 
satisfies for 0 < t < T and a(t) < T - , 

where a(t) = 2t(l - for all t £ [0, T], 

t ' s 

Li(t,x 0 ,y 0 ) = J Lfi(s,x 0 ,y 0 , J G^rXg^Xcyo) dr, 

0 — co 

i>(s) r 
, J ^(T^cy^dT) --|/i(5,x 0 ,y 0 , 

a(s) 0 

, J G 1 (5,T)£f 1 (T,x 0 ,y 0 )dT, J g 1 {T,XQ l yo)dT)ds]ds 

-co a(s) 
t s 

L 2 (t,x 0 ,y 0 ) = j [f 2 (s,x 0 ,y 0 , j G 2 (s,T)g 2 (j,x 0 ,y 0 ) dr, 

0 -co 

b(s) t 

, J £f 2 (T,x 0 ,y 0 )dT) --|/ 2 (5,x 0 ,y 0 , 

a(s) 0 

S 6(5) 

, J G 2 (s,T)g 2 (T,x 0 ,y 0 )dT, j g 2 (T,x 0 ,y 0 )dT)ds]ds 



ds + 



Proof. 

t s 

l|£i(t.*o.yo)ll ^ (1-^) J /i(5,^o.yo. j G 1 (s,T)g 1 (j,x 0 ,y 0 )dT, 

0 -co 

6(s) 

, j g 1 {T l XQ l y Q )dT) 

r s b 

+ fi(s,x 0 ,y 0 , j G 1 (s,T)g 1 (r,x 0 ,y 0 )dT, j g^.x^y^dr) 

t -co a(s) 

t T 

<(1-^ j M.ds + ^j M,ds 
o t 

^[(l-^t + ^r-t)] 

so that: 

Ui(t,x 0 ,y 0 )\\<M ia (t) -(1.13) 
And similarly, we get also 

U 2 (t,x 0 ,y 0 )\\<M 2 a(t) -(1.14) 
From ( 1.13) and ( 1.14), then the inequality ( 1.12) satisfies for all 

0 < t < T and a(t) < T - . m 

II. Approximation Solution Of (Ii) And (I 2 ) 

In this section, we study the approximate periodic solution of (Ij) and 
(I 2 ) by proving the following theorem: 

Theorem 2.1. If the system (10 and (I 2 ) satisfies the inequalities ( 1.2), ( 1.3), 1.4), (1.5),(1.6) and 
conditions(1.7),(1.8) has periodic solutions x = x{t, x 0 ,y 0 ) and y = y(t,x 0 ,y 0 ), then the sequence of 
functions: 

t 

* m +i(t,*o,yo) = ^o(t) + j [fi(s,x m (s,x 0 ,y 0 '),y m (s,x 0 ,y 0 '), 
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, j G 1 (5,T)£f 1 (T,x m (T,x 0 ,y 0 ),y m (T,x 0 ,y 0 ))dT, 

— oo 

i(s) T 

, j g 1 (T,x m (T,x 0 ,y 0 ),y m (T,x 0 ,y 0 ))dT)--jf 1 (s,x m (s,x 0 ,y 0 ), 

a(s) 0 

s 

,ym(.s,x 0 ,y 0 ), j G 1 (s,T)g 1 (T,x m (T,x 0 ,y 0 ),y m (T,x 0 ,y 0 ))dT, 

-co 

b(.s) 

, j £f 1 (T,x m (T,x 0 ,y 0 ),y m (T,x 0 ,y 0 ))dT)ds]d5 •••(2.1) 



a (s) 

with 

XoO-Xo.yo) = F o(0 = *o . for all m = 0, 1,2, • 



ym+i(t,x 0 ,y 0 ) = G 0 (t) + J [Ms,x m (s,x 0 ,y 0 ).y m (s,Xo.yo). 

0 

s 

, j G 2 (s,T)g 2 (j,x m (T,x 0 ,y 0 ),y m (T,x 0 ,y 0 ))dT, 

—00 

, J £f 2 (T,x m (T,x 0 ,y 0 ),y m (T,x 0 ,yo))dT) -- j f 2 (s,x m (s,x 0 ,y 0 ), 

a (5) 0 

s 

,y m (s,x 0 ,y 0 ), J G 2 (s,T)£f 2 (T,x m (T,x 0 ,y 0 ),y m (T,x 0 ,y 0 ))dT, 

-co 

, j g2(T,x m (j,x 0 ,y 0 ),y m (T,x 0 ,y 0 ))dT)ds]ds •••(2.2) 

a (5) 

with 

yo(t.*o>yo) = G 0 (t) = y 0 , for all m = 0,1,2, - 

periodic in t of period T, and uniformly converges asm-) 00 in the domain: 

(t, x 0 , y 0 ) e [0, T] X Dfo x% -(2.3) 

to the limit functions x°(t, x 0 ,y 0 ), y°(t, Xo<yo) defined in the domain ( 2.3) which are periodic in t of period T 

and satisfying the system of integral equations: 
t 

x(f,x 0 ,y 0 ) = F 0 (t) + j [A(s,x(s,Xo,yo),y(s,x 0 ,yo), 
0 

s 

, j G 1 {s,x)g 1 (x,x{x,x 0 ,y 0 ),y{x,x 0 ,y 0 ))dx, 

—00 

b(s) T 

, j g 1 (r,x<j,x 0 ,y 0 ),y(T,x 0 ,y 0 ))dT)--jf 1 (is,x<is,x 0 ,y 0 ), 

a (5) 0 

s 

.y(s,x 0 ,y 0 ), j G 1 (s,T)g 1 (T,x(T,x 0 ,y 0 ),y<j,x 0 ,y 0 ))dT, 

-co 

, j g 1 (T,x(T,x 0 ,y 0 ),y(T,x 0 ,y 0 ))dT)ds]ds •••(2.4) 



and 

t 

y(t.^o.yo) = G 0 (t) + J [/ 2 (5,x(5,x 0 ,y 0 ),y(5,x 0 ,yo), 
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, j G 2 (s,T)g 2 (T,x(T,x 0 ,y 0 ),y(j,x 0 ,y 0 ))dT, 

— oo 

i(s) T 

, j g 2 (T,x(j,x 0 ,y 0 ),y(T,x 0 ,y 0 ))dT) -- j f 2 (s,x(s,x 0 ,y 0 ), 

a (s) 0 

s 

,y(s,x 0 ,y 0 ), j G2(5,T)£f 2 (T,x(T,x 0 ,yo).y(T,^o.yo))dT, 

— CO 

./*(/) 52(T,x(T,x 0 ,y 0 ),y(T,x 0 ,y 0 ))dT)ds]ds -(2.5) 



provided that: 

,0 



(\\x°(t,x 0 ,y 0 )-x m (t,x 0 ,y 0 )\\\ m - 

l||y°(t.^yo)-y I n(^o.yo)llJ- <30 (£_(?o) C ° "' U6) 



for all m > 0 , where 

C 0 = r , £ is identity matrix. 
\ 2 M 2 ) 

Proof. Consider the sequence of functions *i(t, x 0 ,y 0 ),x 2 (t, x 0 ,y 0 ), 
x m (t,x 0 ,y 0 ),-" and y x (t, x 0 ,y 0 ),y 2 (t, x 0 ,y 0 ), ••• ,y m (t, x 0 ,y 0 ), " defined 

by the recurrences relations ( 2.1) and ( 2.2). Each of these functions of sequence defined and continuous in the 
domain (1.1) and periodic in t of period T. 

Now, by the lemma 1.1, and using the sequence of functions (2.1), 
when m — 0, we get: 

By mathematical induction and lemmal.l , we have the following inequality: 

\\x m (t,x 0 ,y 0 )-x 0 \\ < T - M t -(2.13) 

i.e. x m (t,x 0 ,y 0 ) £ £),forall t £ R 1 ,x 0 £ D fl ,y 0 £ D lfv m = 0,1,2, 
and 

T 

\\y m (.t,x 0 ,y 0 )-y 0 \\ < - M 2 -(2.14) 

i.e. y m (t,x 0 ,y 0 ) £ D, for all t £ R 1 ,x 0 £ £>/i,y 0 e D 1/2 ,m = 0,1,2,-, 
Now, from ( 2.13), we get: 

T Y 

\\z m (.t,x 0 ,y 0 )-z Q (t,x 0 ,y 0 )\\ < -(-) (R^ + R 2 M 2 ) -(2.15) 

and 

r 

l|w m (t,x 0 ,yo) - w 0 (t,x 0 ,yo)ll ^ ^HR^ + R 2 M 2 ) - ( 2.16) 

for all t £ R\x 0 £ f/^yo G Di/ 2 , z 0 (t,x 0 ,y 0 ) £ Z^/i^d 
w 0 (t,x 0 ,y 0 ) £ D 3/l , 

i.e. z m (t,x 0 ,y 0 ) £ D 2 and w m (t,x 0 ,y 0 ) £ D 3 , forallx 0 £ f/^yo G £>i/ 2 , 
where 

t 



i(t,^o-yo)= j G 1 (t,s)g 1 (s,x m (s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ))ds 



and 

w m (t,x 0 ,y 0 ) = j £f 1 (5,x m (5,x 0 ,y 0 ),y m (5,x 0 ,y 0 ))d5, m = 0,1,2,-, 

a(t) 

Also from (4.2.14), we get: 

7 8 

\\u m (t,x 0 ,y 0 ) - u 0 (t,x 0 ,y 0 )\\ < -(— ) (H^ + tf 2 M 2 ) - ( 2.17) 

Z A 2 

and 

\\V m {t,x 0 , yo ) -V 0 {t,x 0 ,y 0 )\\ < T -h{H 1 M 1 + H 2 M 2 ) -(2.18) 
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for all t £ R x ,x 0 £ D fl ,y 0 £ D lf2 , u 0 (t,x 0 ,y 0 ) £ D 2 / 2 and 
v 0 (t,x 0 ,y 0 ) £ % 2 , 

i.e.u m (t,x 0 ,y 0 ) £ £> 2 and v m (t,x 0 ,y 0 ) £ £> 3 ,forallx 0 £ D fl ,y 0 £ D lf2 , 
where 

t 

u m (t,x 0 ,y 0 ) = J G 2 (t,5)£f 2 (5,x m (5,x 0 ,y 0 ),y m (5,x 0 ,y 0 ))d5 



and 



6(0 



a(0 



m = 0,1,2,"-, 



Next, we claim that two sequences{x m (t,x 0 ,yo)}m=o<{ym(t<*o<yo)}m=o 
are convergent uniformly to the limit functions x,y on the domain ( 2.3). 
By mathematical induction and lemma 1.1, we find that: 



\\x m +i(t,x 0 ,y 0 ) - x m (t,x 0 ,y 0 )\\ < 



.Y 



< a(t)[^i +Ri(t~ K 3 + hK 4 )]\\x m (t,x 0 ,y 0 ) - x m _ 1 (t,x 0 ,y 0 )\\ 
A \ 

+a(t)[K 2 + K 2 (f K 3 + hK 4 )]\\y m (_t,x 0 ,y Q ) -y m -i(t,x 0 ,y 0 )\\ 



and 

\\ym+i(t,x 0 ,y 0 ) - y m (t,x 0 ,y 0 )\\ < 



< a(t)[Li + HiC— L 3 + ftL 4 )]||x m (t,x 0 ,y 0 ) - a: m _i (t, x 0 ,y 0 )|| + 



(2.19) 



+a(t)[L 2 + h 2 (— l 3 + w, 4 )]||y m (t,x 0 ,y 0 ) -y m -i(t,^o,yo)ll 

*2 



We can write the inequalities ( 2. 19) and ( 2.20) in a vector form: 
C m +i 

(t,x 0 ,y 0 ) < Q(t)C m (t,x 0 ,y 0 ) 

where 

r ft r -\ = (W x m+i( t ' x o.yo)-x m it,x 0 ,y 0 )\\\ 
m+lL ' °' y ° ; Mly m+ i(t.*o.yo)-y m (t.*o.yo)llJ' 

and 

/ o(t)[Ki + ffi(f tf 3 + WQ)] a(t)[/f 2 + « 2 (f ^ + ^ 4 )] N 

<K0 = ^ / 

ya(t)[Li +//i(^L 3 + W. 4 )] «(t)[L 2 + #2(^3 + 

r r >, /llx m (t.*o.yo) -*m-i(t.*o.yo)ll\ 

mL ' °' y ° j V||y m (t,x 0 ,y 0 )-y m _ 1 (t,x 0 ,y 0 )||; 

Now, we take the maximum value to both sides of the inequality ( 2.21), 

for all 0 < t < T and a(t) < T - , we get: 
where Q 0 = max t€[0J] \Q(t)\ . 



Qo = 



T [L, + H^L, + hi,)] T -[L 2 + // 2 (£i 3 + hi,)} 



2 L 1 1V A 2 3 2 

By iterating the inequality (4.2.22), we find that: 

Cm+i — Qo m Co > 
which leads to the estimate: 



Since the matrix Q 0 has maximum eigen-values of (4.1.13) and the series 
(4.2.24) is uniformly convergent, i. e. 



(2.21) 



(2.22) 



■ (2.23) 
(2.24) 



(2.20) 
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lim Y Qj- 1 C 0 = Y CV" 1 C 0 = (E- QoT 1 ^ 

m-»co / i / i 

i=l i=l 

The limiting relation (4.2.25) signifies a uniform convergence of the sequence 
{x rn (t,x 0 ,y 0 '),y m (t,x 0 ,y 0 ')}'^ =0 in the domain (4.2.3) asm^co. 
Let 

lim x m (t,x 0 ,y 0 ) = x°(t,x 0 ,y 0 ) , 

m-»co 

lim y m (t,x 0 ,y 0 ) = y°(t,x 0 ,y 0 ) . 

m-»co / 

Finally, we show that x(t,x 0 ,y 0 ) = x°(t,x 0 ,y 0 ) ED and 

y(f,x 0 ,y 0 ) =y°(t,x 0 ,y 0 ) EDi.forall x 0 £ D fl andy 0 G D lfl . 

By using inequalities (2.1) and ( 2.4) and lemma 1.1, such that: 
t 

j [fi 0> x m (s, x 0 ,y 0 ),y m (s, x 0 ,y 0 ),z m (s, x 0 , y 0 ), w m (5, x 0 , y 0 )) - 
0 

r 

- ^ J A (5- *m (5- *o> yo). y m 0> x 0 ,y 0 ), z m (s, x„, y 0 ), w m (5, x 0 , y 0 ))ds]ds - 
0 

t 

- J [AO, x(s, x 0 ,y 0 ),y(s, x 0 ,y 0 ), z{s, x 0 ,y 0 ), w(s, x 0 ,y 0 )) - 
0 

T 

-j j AO>*0> x 0 ,y 0 ),y(s, x 0 ,y 0 ),z(_s, x 0 ,y 0 ),w(s,x 0 ,y 0 ))ds]ds 



■(2.25) 



(2.26) 



< a(t) (J/fi + ^1(7^3 + ^4)]lkm(s.x 0 ,y 0 ) -x(_s,x 0 ,y 0 )\\ + 

+ [K 2 +R 2 (.j~ i K 3 + hK 4 )]\\y m (s,x 0 ,y 0 ) - y(s,x 0 ,y 0 )\\) < 
T I V 

< -([^ + R^-^ + hK 4 )]\\x m (t,x 0 ,y 0 ) - x(t,x 0 ,y 0 )\\ + 

+ [K 2 + ^2(^3 +hK 4 )]\\y m (t,x 0 ,y 0 )-y(t,x 0 ,y 0 )\\) . 
From inequality ( 2.26) and on the other hand suppose that:. 
\\x m (t,x 0 ,y 0 ) - x(t,x 0 ,y 0 )\\ < e t , 
\\y m (t, x 0 , y 0 ) -y(t,x 0 , y 0 ) || < e 2 . 

Thus 

t 

j [fds,x m (s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ),z m (s,x 0 ,y 0 ),w m (s,x 0 ,y 0 )) - 
0 

T 

-fj fds,x m (s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ),z m (s,x 0 ,y 0 ),w m (s,x 0 ,y 0 ))ds]ds - 
0 

t 

- J [AO, x{s, x 0 ,y 0 ),y(s, x 0l y a ), z(s, x 0 ,y 0 ), w(s, x 0 ,y 0 )) - 
0 

T 

-f j AO>*0> ^o,yo),y0,^o,yo),z0,^o,yo),vi/(s,xo,yo))d5]ds 

0 

- \ [Kl + Rl( ^ 3 + hK4)]€l + 1 2 [K2+ ^T^ 3 + hKi)]€2 

Putting e l — j y and e 2 — j y and substituting in the last equation, we have: 



[K 2 +R 2 (.jf 3 +hK i )] 



j [f t (5, x m (5, x 0 , y 0 ), y m (s, x 0 , y 0 ) , z m (s, x 0 , y 0 ), w m (s, x 0 , y 0 )) - 
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i 

-jj fds,x m (s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ),z m (s,x 0 ,y 0 ),w m (s,x 0 ,y 0 ))ds]ds - 
o 

t 

- j LAO. *0. x 0 ,y 0 ), y(s, x 0 ,y 0 ), z(s, x 0 ,y 0 ), w(s, x 0 , y 0 )) - 
o 

T 

-j j AO.*0. x 0 ,y 0 ),y(s,x 0 ,y 0 ),z(s,x 0 ,y 0 ),w(s,x 0 ,y 0 ))ds]d^ 
o 

< I [Ki + «i(f *3 + WQ)] T 5 + 

z Al 7^i+fii(^ 3 + WQ)] 

+ 7 [K 2 + RiCj- K 3 + hK 4 )] = ^ 

< <?3 + e 4 , 
and choosing e 3 + e 4 — e , we get: 
t 

J LA 0. *m 0. x 0 ,y 0 ), y m (s, x 0 ,y 0 ), z m (s, x 0 ,y 0 ), w m (s, x 0 ,y 0 )) - 
o 

T 

-fj fi(s,x m (_s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ),z m (s,x 0 ,y 0 ),w m (_s,x 0 ,y 0 ))ds]ds - 
o 

t 

- j [A (5, x(s, x 0 ,y 0 ), y(s, x 0 ,y 0 ), z(s, x 0 ,y 0 ), w(s, x 0 , y 0 )) - 
o 

r 

j fi(s,x(s,x 0 ,y 0 ),y(s,x 0 ,y 0 ),z(s,x 0 ,y 0 ),w(s,x 0 ,y 0 ))ds]ds 
o 

for all m > 0, 

t 

i.e. lim [A(5,x m (5,x 0 ,y 0 ),y m (s,x 0 ,y 0 ),z m (5,x 0 ,y 0 ), w m (s,x 0 ,y 0 )) - 

m->co J 
0 

r 

J" /i(5,x m (5,Xo,yo),y m (5,Xo,yo),z m (s,Xo,yo),w m (s,x 0 ,yo))ds]ds = 

0 

t 

- 1 LAO. *0. ^o.yo). yO. *o.yo). z O. *o.yo). w O. *o. yo)) - 

0 

T 

~ffh 0. *0. *o. yo). y 0. *o. yo). ^0. x 0 , y 0 ), w(s, x 0 . y 0 ))ds]ds. 

0 

So x(t, x 0 ,yo) £ £>, and x(t,x 0 ,y 0 ) = x°(t, x 0 ,y 0 ) is a periodic solution of ( II), 

Also by using the same method above we can prove y(t,x 0 ,y 0 ) £ D x , and y(t,x 0 ,y 0 ) = y°(t,x 0 ,y 0 ) is 
also periodic solution of (7 2 ). 

Theorem 2.2. With the hypotheses and all conditions of the theorem 2.1, the periodic solution of integral 

equations ( /i ) and ( I 2 ) are a unique on the domain (1.3). 

Proof. Suppose that x(t, x 0 , y 0 ) and y (t, x 0 , y 0 ) be another periodic solutions 

for the systems (4) and ( I 2 ) defined and continuous and periodic in t of period T, this means that: 
t 

x(.t,x 0 ,y 0 ) = F 0 (t) + J [A0.*0.*o.yo).y0.*o.yo). 

o 

s 

, j G^s.^g^T, x(t, x 0 ,y 0 ),y(r,x 0 ,y 0 )) dT, 



< e 
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bis) 



ats) 

and 



, j g 1 (T,x(T,x 0 ,y 0 ),y(r,x 0 ,y 0 ))dT)-^;jf 1 (s,x(s,x 0 ,y 0 ), 

a(s) 0 

s 

.yO,x 0 ,y 0 ), j G i(. s ' T )9i( T 'X(j,x 0 ,y 0 ),y(T,x 0 ,y 0 ))dT, 

-co 

b(s) 

, j £f 1 (T,f(T,x 0 ,y 0 ),y(T,x 0 ,y 0 ))dT)d5]ds •••(2.27) 



y(t,x 0 ,y 0 ) = G 0 (t) + j [/ 2 (s,x(5,x 0 ,y 0 ),y(s,x 0 ,y 0 ), 
o 

s 

, j G 2 (s,T)g 2 (T,x(T,x 0 ,y 0 ),y(j,x 0 ,y 0 ))dT, 

— oo 

i>(s) t 

, j g 2 (j,x(j,x 0 ,y 0 ),y(j,x 0 ,y 0 '))dT) - - j f 2 (s,x(s,x 0 ,y 0 '), 

a(s) 0 

s 

,9(s,x 0 ,y 0 ), j G 2 {s,T)g 2 {T l x{T,XQ,yQ) l y{T,XQ l yQy)dT, 



bis) 



, j £f 2 (T,f(T,x 0 ,y 0 ),y(T,x 0 ,y 0 ))dT)ds]ds •••(2.28) 



a(s) 

For their difference, we should obtain the inequality: 

\\x(t,x 0 ,y 0 ) -x(t,x 0 ,y 0 )\\ < 

T Y 

< 2 [^i + + Wf 4 )]lk(t.«o.yo) - x(t,x 0 ,y 0 )\\ + 

+ \[K 2 + R 2(£^3 + wr 4 )]||y(t.*b.y 0 ) -Kt.«b.y 0 )ll -( 2-29) 

And also 

\\y(s,x 0 ,y 0 ) - y(s,x 0 ,y 0 )\\ < 

t 

<(1-^)J [Lil|x(s,x 0 ,y 0 )-x(5,x 0 ,y 0 )|| + L 2 ||y(5,x 0 ,y 0 ) - y(5,x 0 ,y 0 )|| 
o 

5 

+ ^-L 3 (H 1 \\x(s,x 0 ,y 0 ) - x(s,x 0 ,y 0 )\\ + H 2 \\y(s,x 0 ,y 0 ) - y(s,x 0 ,y 0 )\\) + 
+hL 4 (H 1 \\x(s,x 0 ,y 0 ) - x(s,x 0 ,y 0 )\\ + H 2 \\y(s,x 0 ,y 0 ) - y(s,x 0 ,y 0 )\\)]ds 

T 

+ Y j"[Lil|x(s,x 0 ,y 0 ) - x(5,x 0 ,y 0 )|| + L 2 \\y(s,x 0 ,y 0 ) - y(s,x 0 ,y 0 )\\) + 
t 

5 

+ ^-L 3 (H 1 \\x(s,x 0 ,y 0 ) -x(s,x 0 ,y 0 )\\ + H 2 \\y(s,x 0 ,y 0 ) - y(s,x 0 ,y 0 )\\) + 

A 2 

+hL 4 (H 1 \\x(s,x 0 ,y 0 ) - x(s,x 0 ,y 0 )\\ + H 2 \\y(s,x 0 ,y 0 ) - y(s,x 0 ,y 0 )\\)]ds 
5 

< a(t) [Li +H 1 {— L 3 +hL i )]\\x(t,x 0 ,y 0 ) -x(t,x 0 ,y 0 ) \\ + 

A 2 

+a{t)[L 2 + H 2 (^L 3 + hL 4 )]\\y(t,x 0 ,y 0 ) - y(t,x 0 ,y 0 )\\ 
so that: 

\\y(t,x 0 ,y 0 ) -y(t,x 0 ,y 0 )\\ < 

T 8 

< -\L X +//!(— L 3 + hL 4 )]||x(t,x 0 ,y 0 ) -x(t,x 0 ,y 0 )\\ + 

+ T -[L 2 + #2(^3 + hL 4 )]\\y(.t,x 0 ,y 0 ) - y(.t,x 0 ,y 0 )\\ - ( 2.30) 
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and the inequalities (4.2.29) and (4.2.30) would lead to the estimate: 
f\\x(t,x 0 ,y 0 ) - x(t,x 0 ,y 0 )\\\ < t\\x(t,x 0 ,y 0 ) - x(t,x 0 ,y 0 )\\\ 
\\\y(t,x 0 ,y 0 ) - y(t,x 0 ,y 0 )\\) ~ 0 \\\y(t,x 0 ,y 0 ) - y(t,x 0 ,y 0 )\\) 

•••(2.31) 

By iterating the inequality (4.2.27), which should find: 
(\\x(t,x 0 ,y 0 ) -x(t,x 0 ,y 0 )\\\ < m (\\x(t,x 0 ,y 0 ) - x(t,x 0 ,y 0 )\\\ 
V||y(t,x 0 ,y 0 ) -y(t,x 0 ,y 0 )\\) ~ 0 V||y(t, x 0 , y 0 ) - y(t,x 0 ,y 0 ) \\) ' 

But Q 0 m -» 0 as m -» oo, so that proceeding in the last inequality which is contradict the supposition It follows 
immediately x(t,x 0 ,y 0 ) = x(t,x 0 ,y 0 ) and 
y(t,x 0 ,y 0 ) = y(t,x 0 ,y 0 ). 

III. Existence of Solution of (1^ and (/ 2 ) 

The problem of existence of periodic solution of period T of the system 
( Ii) and ( l 2 ) are uniquely connected with the existence of zero of the functions A 1 (0,x 0 ,y 0 ) = Ax and 
A 2 (0,x 0 ,y 0 ) = A 2 which has the form: 
A 1 :D fl xD lf2 -^R n 

T 

M0,x 0 ,y 0 )= ^| f 1 (t,x 0 (t,x 0 ,y 0 ),y°(t,x 0 ,y 0 ), 

-co 

bit) 

, j g 1 (s,x 0 (s,x 0 ,y 0 ),y°[[s,x 0 ,y 0 '))ds)dt •••(3.1) 



adt) 

A 2 :D fl xD lf2 

T 

A 2 (0,x 0 ,y 0 ) = -j f 2 (t,x 0 (t,x 0 ,y 0 ),y°(t,x 0 ,y 0 ), 
o 

t 

, J G 2 (t,s)£f 2 (s,x 0 (5,x 0 ,yo),y 0 (5,x 0 ,yo))d5 

-co 

bit) 



, / 92 (,A,*„.y„),A,^„))<*)<« -.(,2, 



a(t) 

where the function x°(t, x 0 ,y 0 ) is the limit of the sequence of the functions x m {t, x 0 ,y 0 ) and the function 

y° (t, x 0 , y 0 ) is the limit of the sequence of the functions y m (t, x 0 , y 0 ). 

Since this two functions are approximately determined from the sequences of functions: 



A lm :D h xD lf2 ^R n 

T 

A lm (0,x 0 ,y 0 ) = - j A(t,x m (t,x 0 ,y 0 ),y m (t,x 0 ,y 0 ), 



L 

, j G 1 (t,s)g 1 (s,x m (s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ))ds 

— CO 

bit) 

, j g 1 (s,x m (s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ))ds)dt •••(3.3) 



ait) 

A 2m :D h xD lf2 ^ R n 

T 

A 2m (0,x 0 .yo) = j;j f2(t,x m (t,x 0 ,y 0 ),y m (t,x 0 ,y 0 ), 

0 

t 

, j G 2 (t,s)g 2 (s,x m (s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ))ds 



www. ij eij ournal. com 



Page I 10 



Periodic Solutions for Non-Linear Systems of Integral Equations 



b(f) 



, j g 2 (s,x m (s,x 0 ,y 0 ),y m (s,x 0 ,y 0 ))ds)dt •••(3.4) 

a(t) 



for all m — 0,1,2, ■ 



Theorem 3.1. If the hypotheses and all conditions of the theorem 2.1 and 2.2 are satisfied, then the following 
inequality satisfied: 

l|Ai(0,x 0 ,y 0 ) - A lm (0,x 0 ,y 0 )|| < d m ••• (3.5) 

l|A 2 (0,x 0 ,y 0 )-A 2m (0,x 0 ,y 0 )|| < r\ m -(3.6) 
satisfied for all m > 0 , x 0 £ D fl and y 0 £ D lf2 , 
where 

d m = ({[K, + R,(^K 3 + hK 4 )] [K 2 + R 2 (^K 3 + Wf 4 )]), 
.Qo^dE-Qor'Co) 

and 

n m = (([Li + HiC^La + W. 4 )] [L 2 + // 2 (^L 3 + «. 4 )]) 

.Qo^i-Qo)- 1 ^). 
Proof. By using the relation ( 3.1) and ( 3.3), we have: 

l|Ai(0,x 0 ,yo)-A 1?n (0,Xo,yo)ll 

< [Ki +«i(f ^3 +hK 4 )]\\x°(t,x 0 ,y 0 ) -x m (t,x 0 ,y 0 )\\ + 

+ [K 2 + R 2 (f ^ 3 + hK 4 )]\\y°(t,x 0 ,y 0 ) - y m (t,x 0 ,y 0 )\\ 

< <([*! + ^3 + Wf 4 )] [*2 + R 2 (f ^3 + WQ)]) 

,(? 0 m+l (£-(?or 1 Co> = d m 

And also by using the relation ( 3.2) and ( 3.4), we get: 
l|A 2 (0,x 0 ,yo)-A 2m (0,Xo,yo)|| - 

< ({[L, + H.C-L, + hL 4 )] [L 2 + H 2 (-L 3 + hL 4 )]j, 

,Q 0 m+ \E-Q o y 1 C 0 ) = Vm 
where (. ) denotes the ordinary scalar product in the space R n . m 



Theorem 3.2. Let the vector functions f x [t, x, y, z, w), f 2 (t, x, y, u, v), 
g t (t, x, y)and g 2 (t, x, y) be defined on the domain: 
G = {0 < 5 < t < T,a < x,y < b,c < z,u < d,e < w,v < /} £ R 1 , 
and periodic in t of period T. 

Assume that the sequence of functions (4.3.3) and (4.3.4) satisfies the 
inequalities: 

T min T A lm (0,x 0 ,y 0 ) <-d m , \ 

a + jMi <x o <b — 2"Mj 

T T 
c+2-M 2 <yo<d-2 M 2 

T max t A lm (0,x 0 ,y 0 ) > d m , 

a + Tj-Mi <x o ^> - jMi 

T T 
c+jM 2 <yo<si- 2 M2 

T min T A 2m (0,x 0 ,y 0 ) < -r| n 

a + Tj-Atfi <x o ^> - jMi 
T T 

c+2M 2 <yo<d- 2 M2 

T max t A 2m (0,x 0 ,y 0 ) >il m 

a + 2"Mj <x o <b — 2"Mj 

T T 
c+2-M 2 <yo<d-2 M 2 



for all m > 0, where 

d m = ({[K, + R.C^-K, + hK 4 )] 



[K 2 +R 2 (yK 3 + hK 4 )]y 



(3.7) 



( 3.8) 



and 
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n m = <([Li + tf^Ls + hL ^ l L 2 + #2(^3 + W. 4 )l) 



,<?o m2+1 (^-<?o)- 1 Co> 
Then the system (4.1.1) and (4.1.2) has periodic solution of period T 

x = x(t,x 0 ,y 0 ) andy = y(t,x 0 ,y 0 ) for which x 0 £ [a + T -M t , b - T -M t ] and 

T T 
y 0 £ [c + -M 2 ,d--M 2 ]. 

T T 

Proof. Let x lt x 2 be any two points in the interval [a + -M 1 ,b — -M^ and y x , y 2 be any two points in the 

T ,. , T 



interval [c + -M 2 , d - ~M 2 ], 
such that: 



AimCd*!^) = T min t A lm (0,x 0 ,y 0 ) 

a + jM\ <x 0 ^> - j M l 

T T 
c+jM 2 <yo<A- 2 M2 

Ai m (0,X2,y 2 ) = T max t A lm (0,x 0 ,y 0 ) 

a+2-Mi<xo^>-2 M l 

r r 

c + jM 2 <yo<d- 2 M 2 

A 2 m(0,x 0 ,y 0 ) , 
^2m(0,x 0 ,y 0 ) , 



( 3.9) 



(3.10) 



(3.11) 



A 2m (0^i,yi)= T min t 

a+^-MjSxoSb-^l 

r r 

c+2-M2<yo<d-2 M 2 

A 2m (0,x 2 ,y 2 ) = r max r 

a+2-Atf 1 <xo^>-2* f l 

r r 

c+2-M2Syo<d-2 M 2 

By using the inequalities ( 3.5), ( 3.6), ( 3.7) and (3.8), we have: 
A 1 (0,x 1 ,y 1 ) = A lm (0,x 1 ,y 1 ) + [A^O,*^) - A^.x^)] < 0 ,} 
Ai(0,x 2 ,y 2 ) = A lm (0,x 2 ,y 2 ) + [Ai(.0,x 2 ,y 2 ) - A lm (0,x 2 ,y 2 )] <0 J 

A 2 (0,x liyi ) = A 2m (0,x 1 ,y 1 ) + [A 2 (0,x 1 ,y 1 )-A 2m (0,x 1 ,y 1 )] < 0 ,} 
A 2 (0,x 2 ,y 2 ) = A 2m (0,x 2 ,y 2 ) + [A 2 (0,x 2 ,y 2 ) - A 2m (0,x 2 ,y 2 )] < 0 J 

•••(3.12) 

It follows from the inequalities (3.11) and ( 3.12) in virtue of the 
continuity of the functions A x (0, x 0 , y 0 ) and A x (0, x 0 , y 0 ) that there exists an isolated singular point (x°, y°) = 
(x 0 ,y 0 ),x° £ [*i,x 2 ] andy 0 £ [yi,y 2 ], 

so that A^O,* 0 ^ 0 ) = 0 and A 2 (0,x°,y°) = 0. This means that the system ( 3.1) and (4.3.2) has a periodic 
solutions x(t,x 0 ,y 0 ),y(t, x 0 ,y 0 ) for which 

x 0 £ [a + ^M^fo-^Mi] andy 0 £ [c + T -M 2 ,d - T -M 2 ]. * 

Remark 3.1. Theorem 3.2 is proved when R n = R 1 , on the other hand as x 0 ,y 0 are a scalar singular point 
which should be isolated (For this remark, see [5]). 

IV. Stability Theorem Of Solution And (/ 2 ) 

In this section, we study theorem on stability of a periodic solution for the integral equations ( 4) and ( I 2 ). 
Theorem 4.1. If the function A 1 (0, x 0 ,y 0 ), A 2 (0,x 0 ,y 0 ) are defined by 

equations ( 3.1) and ( 3.2), where the function x (t, x 0 ,y 0 ) is a limit of the sequence of the functions ( 2.1) , the 
function y°(t,x 0 ,y 0 ) is the limit of the sequence of the functions ( 2.2) , Then the following inequalities yields: 
Ui(0,x 0 ,y 0 )\\<M 1 -(4.1) 
l|A 2 (0,x 0 ,y 0 )|| <M 2 -(4.2) 



and 

||A 2 (0,x 0 1 ,y 0 1 ) - A 2 (0,x 0 2 ,y 0 2 )|| < F^E^F^t) - F 0 2 (t)ll + 

+ ^F 1 F 2 E 2 m + £ 4 ) + ~ 2 EiMGoHt) - G 0 2 (0|| 

•••(4.4) 

for all x°,xl,xl £ D fl , y°,yl,yl £ D lf2 , and a(t) = 2t(l ~\) <\, 
where E 1 = [K x + fl^f ff 3 + hK 4 )],E 2 = [K 2 + R 2 (^-K 3 + hK 4 )], 
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E 3 = [L t + H^L 3 + hLJ], E 4 = [L 2 + H 2 ^L 3 + hL 4 )], 
F 1 = [(1 - - -E,)]- 1 andF 2 = (1 - — N^F.y 1 

Proof. From the properties of the functions x°(t,x 0 ,y 0 ) and y° (t,x 0 ,y 0 ) as in the theorem 4.2.1, the 
functions A 1 = ^(.Xo.yo ), A x = A^Xo.yo ), 

x 0 E f/j^o 6 Dif 2 are continuous and bounded by M 1 ,M 2 in the domain ( 1.3). 
From relation ( 3.1), we find: 

T 

II Ai (0,*o ,y 0 || < -j \\fdt,x\t,x 0 ,y 0 ),y°(t,x 0 ,y 0 ), 
o 

t 

, j G 1 (t,s)g 1 (s,x 0 (s,x 0 ,y 0 ),y°(s,x 0 ,y 0 ))ds, 

-co 

b(t) 

. / ft (,A,W.).y°C S .W.))^)ll* 

a(t) 

by using the Lemma 4.1.1, gives: 
IIAiOUo.yo)!! <M V 

And from relation ( 3.2), we get: 

T 

II a 2 (o,* 0 ,y 0 ll < ^ J ll/2(t^°(t^o-yo).y°(t^o.yo). 

0 

t 

, J G 2 (t,s)5 2 (s,x 0 (s,Xo,yo),y 0 (5,x 0 ,yo))ds, 

-co 

, J £f 2 (s,x 0 (5,Xo,yo),y 0 (5,x 0 ,yo))ds)|| dt 

a{t) 

and using Lemma 4. 1 . 1 , we have: 

l|A 2 (0,Xo,yo)ll <M 2 

By using equation (3.1) and lemma 1.1, we get: 
IIAiCO.xJ.y,, 1 ) - Ad0,xly 0 2 )\\ 

T 

< \ |[/filk°(t,x 0 1 ,y 0 1 ) - x°(t,xly 0 2 )\\ + 

+K 2 \y\t,xly^)-y\t,xly { 2 )\\ + 

+ ^K 3 (R 1 \\x°(t,xly 0 1 ) - x°(t,xly 0 2 )\\ + 

+R 2 \\y\t l xly 0 1 ) - y°(t,x 2 ,y 0 2 )||) + 
+ hK A {R l \\x\t,xly Q 1 )-x\t,xly Q 2 )\\ + 
+R 2 \\y°(t,xly 0 1 )-y°(t,xly 0 2 )\\)]dt 

< [K 1 +R 1 (?-K 3 +hK 4 )] Wx^t.xlyo 1 ) - x°(t,xly 0 2 )\\ + 

A i 

+ [*2 +« 2 (f ^3 +hK,)] lly 0 ^,^ 1 ) -y°(t,x 0 2 ,y 0 2 )|| 

so, 

WA^O.xlyo 1 ) - AdO.xlyo^W < E.Wx^t.xlyo 1 ) - x°(t,xly 0 2 )\\ + 

+E 2 \\y°(t,xly 0 1 ) - y°(t,x 2 ,y 0 2 )|| 

•••(4.5) 

And also by using equation ( 3.2) and lemma 1.1, gives: 
IIAztO^yo^-A.au 2 ,^ 2 )!! 

T 

^^MAt^yo 1 )- At,* 0 2 ,y 0 2 )ll + 

+L 2 \\y\t,xly 0 1 ) - y\t,xi,y 0 2 )\\ + 
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+ -L 3 (H 1 \\x°(t,xly 0 1 )-x°(t,xly 0 2 )\\ + 

+H 2 \\y 0 (t,xly 0 i)-y 0 (t,xly 0 2 )\\) + 
+ hL 4 (H 1 \\x\t,xly 0 1 ) - At,x 0 2 ,y 0 2 )|| + 
+H 2 \\y°(t,x 1 0 ,y 0 1 ) - y°(t,xly 0 2 )\\)]dt 

< [Li + HiC-Ls + ftL 4 )] ||*°(t, *J, V) - x°(t,xly 0 2 )\\ + 
+ ih + H 2 ( ^L 3 + hL 4 )] \\y° (t, xl yo 1 ) - y° (t, x 0 2 , y 0 2 ) || . 



Therefore, 

||A 2 (0,x 0 1 ,y 0 1 )-A 2 (0,x 0 2 ,y 0 2 )|| < E^x^t.xlyo 1 ) - x°(.t,xly 0 2 )\\ + 

+E 4 \\y°(t,xly 0 1 )-y°(t,xly 0 2 )\\ 

•••(4.6) 

where the functions x 0 (t, 4 , y 0 1 ) , y 0 (t, 4 , y 0 1 ), x° (t, 4 , y 0 2 ) and 

y°(t,*o,yo 2 ) are solutions of the equation: 

t 

* 0, 4 . yo ) = *o (0 + f\fi(s,x (s, 4 , y 0 fc ), y (s, 4 , y 0 fc ), 
o 

s 

— CO 

6(5) 

, j g 1 {T,x{x,x^,y^),y{T,x^,y^))dj) - 

T 

-f j fi(s,x(s,x$,y$),y(s, 4.yo). 
o 

s 

, j G^s.T^g^T.xtT.xft.yfo.ytT.xft.yft)) dx, 

-co 

6(5) 

, j g 1 (x,x(x,x^,y^),y(x,x^,y^))dx)ds]ds •••(4.7) 

a (s) 

and 

t 

y(t.4-yo) = G 0 fc (O + | Lf2(5,x(5,4,y 0 fc ),y(5,4,y 0 ; ), 

0 

s 

, j G 2 (s,x)g 2 (x,x(x,x$,y$),y(x,x$,y$)) dx, 

-co 

6(s) 

, j g 2 (j,x(T,x$,y$),y(T,x$,y$))dT) - 

a (5) 

r 

J h(s,x(s, 4-yo).y(s.4.yo). 

0 

s 

, j G 2 (s,x)g 2 (x,x(x,x$,yl$),y(x,x$,y$)) dx, 

-co 

6(5) 

, j g 2 (T,x(T,x%,y%),y(j,xft,yfi))dT)ds]ds •••(4.8) 



a (5) 

where fe = 1, 2. 

From (4.4.7), we get: 

Wx^t.xly^-x'it.xly^W 

< 11^(0-^(011 + 
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+o(t)[Ki + Ri(?-K 3 + hK,)] HAt^yo 1 ) - x\t,xly 0 2 )\\ + 
+a(t)[K 2 +R 2 (^K 3 +hK,)] ||y 0 (t,x 0 1 ,y 0 1 )-y°(t,x 0 2 ,y 0 2 )|| 

< 11^(0-^(011 + 

+ 2 [Kl+ Rd i; K3 + kI<4)] "^^'^.yo 1 ) -^°(t,x 0 2 ,y 0 2 )|| + 
+ + ff 2(^s + hK 4 )] \\y\t,xly 0 l ) -y\t,xly 0 2 )\\ 

such that: 

\\x\t,xly^-x\t,xly&\\ < \\F$(t) - F 0 2 (t)|| + 

+ -E 1 \\x 0 (t,x 1 0 ,y 0 1 )-x°(t,xly 0 2 )\\ + 

+ \E 2 \\y\t l xl,y Q 1 )-y\t l xly Q 2 )\\ 

therefore: 

\\x\t,xly 0 i)-x\t,x 2 ) y 0 2 )\\<(l~E 1 yi\\F 0 \t)-F 2 m + 

+ \E 2 {l-\E 1 y 1 \\y\t,xly 0 1 )-y\t,xly 0 2 )\\ 

•••(4.9) 

And also from relation ( 4.8), we have: 

\\y°(t,xly 0 ^-y°(t,xly 0 2 )\\ < 
<||G 0 1 (O-G 0 2 (OII + 

+a(t)[L 1 + //i(-L 3 + hL 4 )] Wx^t.xlyo 1 ) - x°(t,x 0 2 ,y 0 2 )|| + 

+a(t)[L 2 +H 2 (-L 3 + hL 4 )] HyVt.^.yoO-y^t.^.yo 2 )!! < 
<l|G 0 1 (O-G 0 2 (OI| 2 + 

+ +//i( r L 3 +/ii 4 )] HAt^yo 1 ) -x°(t,xly 0 2 )\\ + 

L A 2 

T S 

+ 2 i L 2 + H 2 Cj- L 3 + ftL 4 )] ||y° (t, x 0 \ y, 1 ) " y° (t. * 0 2 , Jo 2 ) II 

so that: 

lly°(t.*3.yo 1 )-y°(t.*g.yo 2 )ll ^ H G o(0 - c 0 2 (OII + 

+ -f 3 ||x 0 (t,x 0 1 ,y 0 1 )-x 0 (t,x 0 2 ,y 0 2 )|| + 

+ ^ 4 ||y 0 (t,x 0 1 ,y 0 1 )-y°(t,x 2 ,y 0 2 )|| 

and hence: 

\\y 0 (t,4,y 0 1 )-y°(t,x$,y 0 2 )\\ < (1 - ^r^CO ~ G 0 2 (OII + 

+\E 3 {i-\E i r 1 \\x\t,xly 0 1 )-x\t,xly 0 2 )\\ 

•••(4.10) 

Now, by substituting inequality ( 4.10) in ( 4.9), we get: 
\\x 0 (t,xly 0 1 )-x°(t,xly 0 2 )\\ < (1 - ^-^(t) - F 0 2 (OII + 

+ |e 2 (1 - ^E 1 )- i [0.-'^Et)- l \\GZ(t) - G 0 2 (t)|| + 

+ ^ 3 (l-^ 4 )- 1 |k 0 (t,x 0 1 ,y 0 1 )-x 0 (t,x 2 ,y 0 2 )||] 

<(l-^ ir 1 llfo 1 (0-F 0 2 (t)|| + 

+ Ie 2 [(1-^E 1 X1 - ^ 4 )]- 1 ||G 0 1 (t) - G 0 2 (t)|| + 

+ ^f 2 f 3 [(i-^ 1 )(i-^4)]- 1 lk 0 (t^o 1 ,yo 1 )-^°(t^o 2 ,yo 2 )ll] 

putting Fi = [(l- T 1 E 1 )(l- T -E 4 )]-\ 



and substituting in the last inequality, we obtain: 
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llAt^yo 1 )- At,* 0 2 ,y 0 2 )ll < (i-^^r'H^CO-^COII + 

+ T -E 2 F 1 \\Gl{t) - G 0 2 (t)ll + + T -^E 2 E 3 F 1 \\x°(t,xly 0 1 ) - x°(t,xly 0 2 )\\] 
as the F 1 (1- T 1 E 4 ) = (1- T 1 E 1 )- 1 

\\x\t,xly 0 1 )-x\t,xly 0 2 )\\<F 1 a-lE 4 )\\F^t)-Fi(t)\\ + 

+ T -E 2 F 1 \\G^(t) - G 0 2 (t)|| + T -^E 2 E 3 F 1 \\x°(t,x 1 0 ,y 0 1 ) - x°(t,xly 0 2 )\\] 

which implies that: 

\\x\t,xly 0 l ) - x\t,xly 0 2 )\\ << F^l - |e 4 )(1 - ^WMl^Ct) - (Oil + 

+ '^E 2 F 1 (1 - ^E^F^WGlit) - G 2 (t) II 

putting F 2 = (1 - — E 2 E 3 F 1 ) 

and substituting in the last inequality, we obtain: 

\\x\t.xiyo 1 ) -x\t,xly 0 2 )\\ < F^l - ^E 4 )\\F^t) - F 0 2 (t)|| + 

+ ^ 2 F 1 F 2 ||G 0 1 (t)-G 0 2 (t)ll -(4.11) 
Also, substituting the inequalities ( 4.1 1) in ( 4.10), we find that: 

Wy^t.xly^-y^it.xly^W < (1 - T -E^\\Glif) - G 0 2 (t)|| + 

+ T -E 3 (l- T 1 E 4 y 1 [F 1 F 2 (l- T 1 E 4 )\\Ft(t)-F 2 (t)\\ + 

+ T -E 2 F 1 F 2 \\Gl{t) - G 0 2 (t) ||] 

and hence 

\\y\t,xly 0 1 )-y°(t,x 2 ,y 0 2 )\\<^E 3 F 1 F 2 \\F^t)-F 2 (t)\\ + 

+ [Fi(l - ^Fi) + ^F!F 2 F 2 ] || Gq 1 (t) - G 0 2 (t) II 

•••(4.12) 

so, substituting inequalities ( 4.1 1) and ( 4.12) in inequality ( 4.5), we get the inequality ( 4.3). 

and the same, substituting inequalities ( 4.1 1) and (4.12) in inequality 
( 4.6), gives the inequality ( 4.4). ■ 
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